Front speed enhancement in cellular flows 
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The problem of front propagation in a stirred medium is addressed in the case of cellular flows 
in three different regimes: slow reaction, fast reaction and geometrical optics limit. It is well known 
that a consequence of stirring is the enhancement of front speed with respect to the non-stirred 
case. By means of numerical simulations and theoretical arguments we describe the behavior of 
front speed as a function of the stirring intensity, U. For slow reaction, the front propagates with a 
speed proportional to (7 1//4 , conversely for fast reaction the front speed is proportional to U 3 ^. In 
the geometrical optics limit, the front speed asymptotically behaves as U/\nU. 
PACS: 



Front propagation in a stirred medium is an 
important problem in a number of fields ranging 
from combustion to plankton dynamics. For a re- 
alistic study of such a class of problems one has to 
take into account the modification of the advect- 
ing flow induced by the reaction, e.g. in combus- 
tion. However, many features can be understood 
by neglecting the back-reaction on the velocity 
field. The problem addressed here is the enhance- 
ment of the front speed induced by a certain class 
of flows. In particular, we consider front propa- 
gation in a two dimensional laminar flow with a 
stationary vortical structure in different regimes, 
namely slow reaction, fast reaction and geometri- 
cal optics limit. This last limit corresponds to a 
very sharp front propagating as an optical front, 
i.e., according to the Huygens principle. We pro- 
vide predictions on the dependence of the front 
speed on the flow intensity, which are confirmed 
by numerical simulations. 



I. INTRODUCTION 

The study of front propagation of a stable phase into 
an unstable one encompasses several issues of great in- 
terest [Q as flame propagation in gases B, population 
dynamics of biological communities (plankton in oceans) 
H and chemical reactions in liquids B . A common fea- 
ture of all these phenomena is that they take place in a 
strongly deformable medium such as a fluid. The inter- 
play among transport, diffusion and reaction is therefore 
a crucial problem with several open issues (e.g. for ques- 
tions concerning combustion see Ref. ||). 

In the most compact model of front propagation the 
state of the system is described by a single scalar field 
0(r,t), that represents the concentration of products. 
The field 9 vanishes in the regions filled with fresh mate- 
rial (the unstable phase), equals unity where only inert 
products are left (the stable phase) and takes intermedi- 
ate values wherever reactants and products coexist, i.e., 



in the region where production takes place. Here we as- 
sume that the concentration of chemicals does not mod- 
ify the underlying flow. Therefore, in the following, we 
consider the velocity field as given. This approximation, 
hardly tenable in the context of flame propagation in 
gases, is rather appropriate for chemical front propaga- 
tion in some liquid solutions • Under these simpli- 
fying assumptions, the evolution of 6 is described by 



d t 9 + u- VO = DA6+-f(0) 

T 



(1) 



where the second term on the l.h.s. accounts for the 
transport by an incompressible velocity field. On the 
r.h.s the first term describes molecular diffusion and 
the second one describes the production process with 
time scale r. We will first consider a production term 
of Fischer-Kolmogorov-Petrovski-Piskunov H (FKPP) 
type, i.e., a function f{9) convex (f"(9) < 0) and posi- 
tive in the interval (0, 1), vanishing at its extremes, and 
f(0) = 1. Here we take f(0) = 0(1 - 0). It is also of 
interest to consider a production term in the form of the 
Arrhenius law, f(0) = (I — 6*)-exp(— c /0), where C is the 
activation concentration. The latter choice is more per- 
tinent to the study of flames and/or chemical reactions 

|0. .... 

Until now we did not specify any details of the veloc- 
ity field. In many engineering applications u is turbulent. 
In this paper we investigate front propagation in laminar 
flows, which, albeit simpler than turbulent ones, show re- 
markable qualitative similarities with more complex flows 
Specifically, we consider a two dimensional station- 
ary incompressible flow with cellular structure (see also 
[fTo| p^| ) u = (—dytp, d x ip) with the streamfunction |L3| 
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We considered L-periodic boundary conditions in y and 
an infinite extent along the x-axis. This kind of flow 
is interesting because, in contrast to shear flows, all the 
streamlines are closed and, therefore, the front propaga- 
tion is determined by the mechanisms of contamination 
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of one cell to the other JlCjQ. Since we are interested 
in the propagation in the ir-direction, the boundary con- 
ditions are set to 9(~oo,y;t) — 1 and 9(+oo,y;t) = 0. 
The maximum principle ensures that at later times the 
field still takes values in the range < 6 < 1 [jlj. The 
instantaneous front speed is defined as 



optics limit. Sect. [v| is devoted to some concluding re- 
marks. The Appendix contains the details of the numer- 
ical method used in the simulations. 



II. UPPER BOUND FOR THE FRONT SPEED 



vt(t) = - 



dxd t 9(x,y;t) . 



(3) 



This expression defines the so-called bulk burning rate 
which coincides with the front speed when the lat- 
ter exists, but it is also a well defined quantity even 
when the front itself is not well defined. The asymp- 
totic (average) front speed, Vf, is determined by Vf = 
liniT^ool/T / dtv { (t). 

In a medium at rest, it is known that Eq. (|l|), for 
FKPP nonlinearity, generates a front propagating, e.g., 
from left to right with an asymptotic speed vq = 2y/ D/t 
and a reaction region of thickness £ = 8 V Dt M. In 
the more interesting case of a moving medium, the front 
will propagate with an average speed Vf greater than vq 
p| , ^2[ . The front velocity Vf is the result of the inter- 
play among the flow characteristics (i.e. intensity U and 
length-scale L), the diffusivity D and the production time 
scale r. The goal of our analysis is to determine the de- 
pendence of Vf on such quantities. In particular, intro- 
ducing the Damkohler number Da = L/(Ut) (the ratio 
of advective to reactive time scales) and the Peclet num- 
ber Pe = UL/D (the ratio of diffusive to advective time 
scales) , we seek for an expression of the front speed as an 
adimcnsional function Vf/v = <j){Da, Pe) > 1. We will 
see that a crucial role in determining such a function is 
played by the renormalization of the diffusion constant 
and chemical time scale induced by the advection Jl4],|l5| . 

Moreover, we consider an important limit case, i.e., 
the so called geometrical optics limit, which is realized 
for {D.t) — > maintaining D/t constant Jl^|. In this 
limit one has a non zero bare front speed, vq, while the 
front thickness £ goes to zero, i.e., the front is sharp. In 
this regime the front dynamics is well described by the 
so-called G-equation J|||[l7| 



dG 
~dt 



u-VG = v \VG\ 



(4) 



The front is defined by a constant level surface of the 
scalar function G(r,i). Physically speaking, this limit 
corresponds to situations in which £ is very small com- 
pared with the other length scales of the problem. Also 
in this case we provide a simple prediction for the front 
speed, which turns out to be expressible as an adimen- 
sional function Vf/vo = tp(U/vo)- 

The paper is organized as follows. In Sect. ^ we dis- 
cuss a theoretical upper bound for the front speed which 
becomes an equality in the limit of (very) slow reaction. 
In Sect. [II we present a numerical study for slow and 
fast reaction, comparing the results with a phenomeno- 
logical model. In Sect. [V we consider the geometrical 



For a generic incompressible flow and a generic produc- 
tion term which has a bounded growth rate, c(0), i.e., 
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sup c(9) = sup — < oo 

ft ft t 



(5) 



it is possible to establish an upper bound for the speed 
of front propagation. Explicitly, we have 



Vf 



< 2^D^< 



(6) 



where D c ff is the effective diffusion coefficient in the re- 
direction, which can be derived from Eq. ([l]) by switching 
off the production term. This bound is the consequence 
of the deeply rooted link existing between front propaga- 
tion and advective transport. 

We start the derivation of Eq. (^) by recalling the fun- 
damental relation among the solution of the PDE (0) and 
the trajectories of particles advected by a velocity field 
u(r, t) and subject to molecular diffusion [pJLpll 



A) = (e{r{Q),Q)e^ c(e(r(s) ' s))ds 



(J) 



The average is performed over the trajectories evolving 
according to the Langevin equation 



dr(t) 
dt 



u(r(i),t) + V2D~v(t) 



(8) 



with ending point r(t) = x. The white noise term 
V2Dr)(t) accounts for molecular diffusion. 
Since the growth rate is bounded, (^) yields the inequal- 
ity 



0(t,x) < (0(r(O),O)), ) exp(c n 



(9) 



In the previous inequality, the term in angular brack- 
ets denotes the probability that the trajectory ending 
at x was initially located at the left of the front in- 
terface. For FKPP production term the maximum oc- 
curs for 8 = 0, i.e., c{6) < c max = c(0) = 1/t. Un- 
der very broad conditions, i.e., nonzero molecular diffu- 
sivity and finite variance of the velocity vector poten- 
tial pfj|-p2[, it is possible to show that asymptotically 
the particles undergo a standard diffusion process with 
an effective diffusion coefficient Z? c ff, always larger than 
the molecular value D. The issue of single particle diffu- 
sion, and the problem of finding the effective diffusiv- 
ity has been the subject matter of a huge amount of 
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work (see e.g. [E2| for a recent review). In the pres- 
ence of an asymptotic standard diffusion, we can sub- 
stitute the term (0(r(O),O)) , with the Gaussian result 

1 - ^erfc(-af/vOT5t) ~ exp[-x 2 /(AD cS t)]/y/2nD cS t 
with exponential accuracy. We thus obtain (9(x, t) < 
exp [c max t — x 2 /(4D e g t)] / \/2T{D C Ht. It is thus clear 
that at the point x the field 9 is exponentially small until 
a time t of the order of x/\/4D c ff c max . We finally obtain 
the upper bound for the front velocity Vf < \/4D c ff c max , 
which is Eq. (ph. 

The bound (O) becomes an equality in the limit of very 
slow reaction. If r is the slowest time scale under con- 
sideration, advection and molecular diffusion act jointly 
to build an effective diffusion process, essentially unaf- 
fected by the reaction. In this case the front width is 
large enough and the reaction takes place in a region 
of effective diffusivity. Therefore, it is allowed to sub- 
stitute Eq. (jl]) with an effective reaction-diffusion equa- 
tion d t 9 = Y,ij D f a ii8+ where D if is thc 
eddy-diffusivity tensor p^] . For the FKPP nonlinear- 
ity, this last equation gives rise to Vf ~ 2-y/D e g/r, where 
D c s — Dff l5] Jl0|Jil[] . One can find a detailed derivation 
of this formula in Ref. |L4| . 

For cellular flows, it is known that D e g ~ \J\JLD 
p3| p5[ . Inserting this expression in Eq. (||) one obtains 
Vf oc C/ 1 / 4 |1C}] , remarkably close to the observed ones for 
Da <C 1 (see next Section). 



III. FRONT SPEED IN THE REACTION 
ADVECTION DIFFUSION EQUATION 



is the rate of change of ©i due to reaction taking place 
within the cell. On the basis of the numerical results (see 
below), we will show that the space-discretized macro- 
scopic reaction-diffusion equation 

-^e, = d cB ( - ©, + le^) + — F(6i) (10) 

dr \2 I J T of f 

is a pretty good model for the front dynamics. The ef- 
fect of the advective field is to renormalize the values 
of the diffusivity, D — > D c g(D, U, L), and the reaction 
time scale, r — > t r(t,U,L). This is why the velocity 
does not appear any more in the effective dynamics, de- 
scribed by Eq. ( [To|) [ ^o|Jll|l . The renormalized diffusivity 
D e ff accounts for the process of diffusion from cell to cell 
as a result of the nontrivial interaction of advection and 
molecular diffusion |23|~|25|]. The renormalized reaction 
time T e ff amounts to the time that it takes for a single 
cell to be filled by inert material, and depends on the 
interaction of advection and production. Of course, also 
the production term will be affected, i.e. / — * F, but, as 
we will see, such a modification is not dramatic. Indeed, 
the modified production term turns always to be in the 
FKPP universality class. 

The limiting speed of thc front in the moving medium 
turns out to be v c ff ~ -Doff/Tcff, similar to Eq.(||) Jl4) . 
The problem is now reduced to derive the expressions for 
the renormalized parameters by means of physical con- 
siderations. 

In the following sections, using as an interpretative 
framework the above described macroscopic model, we 
will present the results of detailed numerical simulations 
for slow (Da <C 1) and fast (Da ^> 1) reaction. 



The bound (|6|) is very general and holds for generic 
incompressible flows and production terms. Here, by 
means of numerical simulations, we consider the front 
propagation problem arising in the reaction advection 
diffusion equation (|lj) for the particular case of the cellu- 
lar flow (|J) of stirring intensity U and FKPP nonlinear- 
ity (with characteristic time, t). In our discussion, we 
always suppose that the diffusion time scale is the slow- 
est occurring one, i.e., L 2 /D <C L/U, r and thus Pe^> 1 
and DaPe < 1. 

Now before presenting the numerical results it is help- 
ful to introduce a macroscopic description of the problem 
which will reduce it to an effective reaction-diffusion pro- 
cess with renormalized coefficients. 

The basic observation is that the dynamics of 9 in 
cellular flows is characterized by the the cell-size L, so 
that we can perform a space discretization that reduces 
each cell, Ci, to a point, i, mapping the domain - a 
two-dimensional infinite strip - onto a one-dimensional 
lattice. The field 6 becomes a function defined on the 
lattice Oi = L~ 2 J c 9dxdy. Integrating Eq. ^ over 

the cell Ci, we obtain = J^+i — Ji + \i where Ji = 
L~ 2 J, ft Dd x 9 dy is the flux of matter through the left 
boundary of the i-th cell, and \i — L~ 2 J c r^ 1 f(9) dxdy 



A. Slow reaction regime 

At small Da, the reaction is significantly slower than 
the advection, and consequently the region where the re- 
action takes place extends over several cells, i.e., the front 
is distributed. To obtain the expression of D Q g we neglect 
the reaction term in Eq. (Q), which reduces to the equa- 
tion for a passive scalar in a cellular flow. This is a well 
studied problem, the solution of which is |^3| : 

^ ~ Pe 1 ' 2 Pe » 1 . (11) 

For large Pe (D small) the cell-to-cell diffusion mech- 
anism can be qualitatively understood in the following 
way. The probability, p, for a particle to jump across the 
boundary of the cell, within a circulation time L/U, by 
virtue of molecular diffusion can be estimated as the ratio 
of the diffusive motion across streamlines, 0(^J DL/U), 
to advective motion along streamlines, O(L). As a re- 
sult p ~ (D /(U L)) 1 / 2 , hence the effective diffusivity 
D cS ~pUL~ DPe 1 ' 2 . 

To obtain the expression of the typical time it takes 
for a whole cell to react, let us concentrate on the re- 
action in a single cell: it is first invaded by a mixture 
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of reactants and products (with a low content of prod- 
ucts, 0; <C 1) on the fast advective time scale; subse- 
quently complete reaction (0j w 1) is achieved on the 
slower time scale r e g ~ r (see Fig. Q). In this regime, 
the front speed is well approximated the homogenization 
result Vf — 2\/D c fi, discussed in the previous section. 




FIG. 1. Six snapshots of the field within the same cell, at 
six successive times with a delay t/6 (from left to right, top 
to bottom), as a result of the numerical integration of equa- 
tion Here Da ~ 0.4, Pe ~ 315. Black stands for = 1, 
white for 8 = 0. 

To check these ideas, we performed numerical simula- 
tions of Eq. (jl]) , with a FKPP production term (see ap- 
pendix for details on the numerical technique) . In Fig. |^, 
we show the result of the calculations for the front speed 
V[ in dependence on Da, the slow reaction corresponds 
to the plateau at Da <C 1. 
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FIG. 2. The ratio of measured front speed, v e g, to the max- 
imal one, v/ 4-D c ff/r, as a function of the Damkohler number, 
Da. For Da <C 1 the front propagates with the maximal 
velocity whereas for Da 3> 1 the speed slows down with in- 
creasing Da. The behavior Da -1 / 2 is drawn for reference. 



B. Fast reaction regime 

We have now to repeat the estimation of D c s and T e s 
for the fast reaction regime, i.e., for large Da. Since we 



work always in the regime of large Peclet numbers, all 
the above arguments for the effective diffusion still hold, 
while the effective time scale is different. At large Da, the 
ratio of time scales reverses, and in a (now short) time 
t two sharply separated phases emerge inside the cell. 
In this regime indeed, the interface is thin compared to 
the cell size. The cell filling process is characterized by 
an inward spiral motion of the outer, stable phase (see 
Fig. ||), at a speed proportional to U, as it usually hap- 
pens for a front in a shear flow at large Da. Therefore, 
the 9 — 1 phase fills the whole cell on the advective time 
scale, giving T c g ~ L/U . 




FIG. 3. Six snapshots of the field 9 within the same cell, at 
six successive times with a delay (L/U)/6 (left to right, top 
to bottom). Here Da = 4, Pe = 315. A spiral wave invades 
the interior of the cell, with a speed comparable to U. 

With respect to the upper bound (^J) we observe for fast 
reaction a significative slowing-down of the front speed 
signaled by a different scaling dependence on the param- 
eters U,D,L,t (see Fig. |). 
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FIG. 4. The renormalized reaction term r/r e gF(B) for 
three different parameter: Da ~ 4 (□), Da ~ 2 (o) and 
Da ~ 0.4 (x). The continuous line shows /(#)• The dotted 
and dash-dotted lines are the slopes (0.2 and 0.4) proportional 
to Da -1 in the region of slow advection. 

We have now to look at the shape of the effective reac- 
tion term t^j F(Q) appearing in the renormalized equa- 
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tion dicf). As shown in Fig. ||, for small Da, the effective 
production term is indistinguishable from the "bare" one. 
Increasing Da, the reaction rate tends to reduce, induc- 
ing the slowing-down of the front speed. For « 0, 
the effective production term essentially coincides with 
the microscopic one. However, there is an intermedi- 
ate regime characterized by a linear dependence on the 
cell-averaged concentration, with a slope proportional to 
Da^ 1 . This is in agreement with a typical effective reac- 
tion time T c s ~ tDcl (see below Eq. ([12|)), To measure 
the macroscopic quantities F(Q) and 0, one simply inte- 
grates numerically both 9 and f(8) at a fixed time over 
a cell volume. 

It is worth to remark that, notwithstanding the change 
of shape of the effective chemical potential, the produc- 
tion term remains in the FKPP universality class. 

Let us stress that what we call fast reaction regime is 
still not in the geometrical optics limit. Indeed, to obtain 
this limit it is not sufficient to take r small, but one has 
to take also D small in such a way that vq ~ 2\JDJt is 
constant when (D, r) goes to zero and the front thickness 
£ oc V Dt is negligible with respect to the cell size. In 
the fast reaction regime studied here the condition on 
the front thickness holds. 



Collecting the information about fast and slow reaction 



results as those of FKPP-type reaction presented above, 
i.e., one has the two scaling laws Vf oc U 1 ' 4 and 
at fast and slow advection respectively (see H). 
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FIG. 5. The front speed Vf as a function of U, the typical 
flow velocity. The lower curve shows data at r = 20.0 (fast 
advection). The upper curve shows data at r = 0.2 (slow 
advection). For comparison, the scalings U 1 ^ and U 3 ^ are 
shown as dotted and dashed lines, respectively. The horizon- 
tal line indicates vo, the front velocity without advection, for 
t = 0.2. 



Tcft 
T 



1 

Da 



Da < 1 
Da > 1 



(12) 



and D c g ~ DPe 1 ^ 2 , we can derive the scaling of the effec- 
tive speed of front propagation for a cellular flow. Indeed, 
recalling that Vf ~ \/D c g/T c g, we have the final result 



V{ ( Pe 1 / 4 

v ~ i Pe^Da- 1 ' 2 



Da < 1 , Pe > 1 , . 

Da» l,Pe> 1 [ ' 



The case of Pe <C 1 is less interesting because the dy- 
namics is dominated by diffusion. 

At small Da the front propagates with an effective ve- 
locity scaling as the upper bound derived above, that is 
as Pe 1 / 4 . At large Da, the front speed enhanced is less 
effective than at small Da: according to Eq. (|l3|), we 
have Vf/vo ~ Da~ x l 2 for Da ^> 1. In terms of the typical 
velocity of the cellular flow, we have Vf oc C/ 1 / 4 for slow 
reaction (U ^> L/t, or equivalently Da -^i 1) whereas 
Vf oc J7 3/4 for fast reaction (U < L/t, or Da > 1). The 
scaling Vf oc C/ 1 / 4 for slow reaction (i.e. fast advection) 
is a consequence of the well known result D c ff oc DPe 1 
p3| in the homogenization limit |l^,Q, it has been ob- 
tained in JlCfl . The numerical results are summarized in 
Fig- 1 

As a remark we mention that, for the class of boundary 
conditions investigated here, where the region of initially 
burnt material extends to infinity, no quenching p6| takes 
place no matter of the used production term. Indeed, 
Arrhenius-type nonlinearity substantially gives the same 



IV. FRONT SPEED IN THE GEOMETRICAL 
OPTICS REGIME 

When the front thickness and the reaction time are 
much smaller than the length and time scales of the ve- 
locity field fluctuations one has the geometrical optics 
regime |l6|,^7j ■ In this case the front is a sharp interface 
separating the reactants from the products, and can be 
modeled in the framework of the G-equation (Q) ||[l6| . 
Physically speaking, one uses the G-equation when the 
front thickness is very thin and it is hard to resolve the 
diffusive scale. 

As far as the cellular flow is concerned, the front bor- 
der is wrinkled by the velocity field during propagation 
and its length increases until pockets of fresh material 
develop |2g|-|3C|] . After this, the front propagates periodi- 
cally in space and time with an average speed Vf, which is 
enhanced with respect to t he p ropagation speed vq of the 
fluid at rest |3l| (see Ref. |p9| for some pictorial views). 

The problem addressed here is the dependence of the 
effective speed Vf on the flow intensity, U, and the bare 
velocity, vq, that is expected of the form 0: 



2 = n>(* 



(14) 



where ip(U) is a function which depends on the flow de- 
tails. 

As far as we know, apart from very simple shear flows 
(for which i)(U) =l+U 016]), there are no methods 
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to compute il>(U) from first principles. Mainly one has 
to resort to numerical simulations and phenomenological 
arguments. 

For turbulent flows, by means of dynamical renormal- 
ization group techniques, Yakhot [B2j proposed 



o(u/v f r 



(15) 



with a = 2. Now U indicates the root mean squared av- 
erage velocity (see also ||[53)). Therefore, from ( |l5| ) one 
has that V{ — > U/yJlii(U) for U — > oo. 

For the cellular flow under investigation, albeit the ex- 
act form of the function ip(U) is not known, a simple 
argument can be given for an upper and a lower bound 
by mapping the front dynamics onto a one-dimensional 
problem. The starting point is the following observa- 
tion. In the optical regime, since the interface is sharp, 
i.e., 0(x,y) is a two-valued function (6 = 1 and 9 = 0), 
we can track the farther edge of the interface between 
product and material (xm(£), 2/m(*))j which is defined 
as the rightmost point (in the x-direction) for which 
0(xm ,VM',t) = 1. Then we can define a velocity 



lim 

t — >oo 



XM{t) 



t 



(16) 



which gives an equivalent value of the standard definition 
within less than 2% (see Fig. m. 




FIG. 6. Front velocity as a function of time, Vf(t) mea- 
sured in the standard way (^) (□), and as in Eq. ( jig ) (o). 
The straight lines represents the average (over a period) of 
the measured value, and the lower and upper bounds. The 
simulation parameters are U = 4, Vo = 1 and L = 2ir. 

In Fig. ^| we show the time evolution of the point 
{%M{t),yM{t))- After a transient, in the unit cell [0, 2vr] 
(we describe the case L = 2ir, i.e., the one adopted here) 
the point 2/m(*)) moves to the right along the 

separatrices of the streamfunction (||), so that yni(t) is 
essentially close to the values or it. Along this path one 
can reduce the edge dynamics to the ld-problem 



da; 



dt 



= v + Uf3\sm(x M (t))\ 



(17) 



where the second term of the r.h.s. is the horizontal 
component of the velocity field. We have neglected the 
y-dependence, replacing it with a constant (3 which takes 
into account the average effect of the vertical component 
of the velocity field along the path followed by (xm, Dm)- 
By solving (|17| ) in the interval xm & (0, 2tt) one obtains 
the time, T, needed for x m to reach the end of the cell. 
The front speed, as the speed of the edge particle, is then 
given by V{ = 2n/T. The final result is 



^ /3 (W)=7rvW) 5 



lln" 



(18) 



Note that (18) is valid only for U(3 > 1. 




FIG. 7. Time evolution of the edge point: xm (t) (solid line) 
and yu{t) (dashed line). The simulation parameters are the 
ones of Fig. ^ 

We have taken j3 = 1 for the upper bound and (3=1/2 
(which is the average of |cos(j/)| between and 7r) for 
the lower bound. We have also computed the average 
of | cos(j/m(£))| in a period of its evolution (see Fig.g) 
obtaining j3 ~ 0.89 which gives indeed a very good ap- 
proximation of the measured curve (see Fig. ||) . We stress 
that the theoretical curve is not a fit, but it just involves 
the measured parameter (3. 

This agreement is an indication that the average of 
|cos(yw(*))| depends on U and vq very weakly (as we 
checked numerically). Previous studies |50| reported an 
essentially linear dependence of the front speed on the 
flow intensity, i.e., Vf oc U for large U which is not too 
far but different from our result. A rigorous bound has 
been obtained in Ref fMfl by using the G-equation: 



Vf > f//(log(l + U/vo)) 



(19) 



As one can see from Fig. g, the lower bound (|19|) seems 
to be closer to the numerical data than the one obtained 
with (3 = 1/2 in (|l|). From Eq. |[|), asymptotically (i.e. 
for U 3> «o) one has: V{ ~ U/ln(U) which corresponds 
to (||) for a = 1. 
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FIG. 8. The measured tp(U/vo) as a function of U/vq (□), 
the Yakhot formula ( |l5| ) with a = 2 (o), the function for 
j3 = 1, 1/2 (dashed and dotted lines) and for /3 = 0.89 (solid 
line). The dashed-dotted line is the bound (|l9|). 

Expressions as (|l5|) have been proposed for flows with 
many scales as, e.g., turbulent flows, and in the literature 
different values of a have been reported |||3j|. The fact 
that also the simple one-scale vortical flow investigated 
here displays such a behavior may be incidental. How- 
ever, we believe that it can be due to physical reasons. 
Indeed, the large scale features of the flow, e.g. the ab- 
sence of open channels (like for the shear flow) can be 
more important than the detailed multi-scale properties 
of the flow @. 

A definite answer to this question is beyond the scope 
of this paper, however it could be an interesting point for 
future investigations. 
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APPENDIX A: NUMERICAL METHOD 

We introduce a lattice of mesh size Ax and Ay (for 
the sake of simplicity we assume Ax = Ay) so that the 
scalar field is defined on the points x n m = (nAx,mAy): 
9 n ,m(t) = 9(nAx,mAy,t). 

Giving the field at time t, the algorithm computes the 
field at time t + At, and the integration method de- 
pends on which kind of physical regime we are interest in: 
reaction-advection-diffusion equation (Q) or optical limit 



1. Reaction Diffusion Equation 

In numerical approaches one is forced to discretize 
the dynamics, so let us consider the case of a velocity 
field which is always zero apart from (5-impulses at times 
t = 0, ±At, ±2At, ±3At, . . . 



We addressed the problem of front speed enhance- 
ment induced by stirring due to a cellular flow in dif- 
ferent regimes. In the slow reaction case a rather gen- 
eral result (based on homogenization techniques) gives 
Vf ~ U 1 / 4 ; for the fast reaction case, physical arguments 
give V{ ~ £/ 3 / 4 . In the geometrical optics limit one finds 
that V{ is a linear function of U, apart from logarithmic 
corrections. All these results has been confirmed by nu- 
merical simulations. 

The steady cellular flow treated here, albeit its sim- 
plicity, provides a paradigm that can be insightful for the 
study of front propagation in more general flows. For in- 
stance, in the geometrical optics limit the asymptotical 
behavior Vf ~ U / ln(?7) is rather similar to the one found 
by Yakhot (v { ~ U/y/ln(U)) fH for turbulent velocity 
fields. 

One could conjecture that there are nontrivial reasons for 
this similarity: the front propagation principally involves 
the largest scales. In this context, it is not surprising that 
a multiscale process (as turbulence) or a single scale pro- 
cess can yield similar results. 



u(x,f) 



71— — OO 



u(x)6(t - nAt) 



(Al) 



In such a case the Lagrangian evolution is given by a 
conservative map (in 2d the map is symplectic due to 
incomprcssibility) 

x(t+ At) = F At (x(i)) . (A2) 

If also the production term is zero apart from ^-impulses 



/(*) 



n— — oo 



g(0)5(t - nAt) 



one can introduce a reaction map 

9(t + At)=G At {9(t)) . 



(A3) 



(A4) 



Let us remark that choosing a 5-impulsed production 
term can be particularly relevant in some experimen- 
tal settings, i.e., when one considers periodic illumina- 
tion in light-sensitive chemical reactions as in Ref. [p7[. 
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The concentration field 0(x, t + At — 0) is obtained from 




or, in other terms: 



(x,t+At) = (G At (e{F 



-A/ 



(x-V2DAt ri(t)),t)} (A6) 



which is equivalent to Eq. ([?]). Let us remark that 
Eq. ( [A6p is exact if both the velocity field and the re- 
action are ^-imp ulse d processes. However one can also 
use the formula (A6) as a practical method for the nu- 
merical integration of Eq. ([l]) if one assumes small enough 
At, so that the Lagrangian and reaction maps are given 
at the lowest order by 



F At (x) ~ x + u(x)A£ , 



FIG. 9. Pictorial scheme of the numerical algorithm for 
RAD systems. Here -q a = V2BAtW° where W Q is a stan- 
dard Gaussian variable. 



To correctly simulate the diffusion process we have to 
impose a relation between D, Ax and At to assu re that 
diffusion transports a particle over distances ~ \J2D At 
much larger than the grid size Ax (see Fig ^ . 



2. Geometrical optics limit 



From an algorithmic point of view the whole process 
between t and t + At, Eq. (|A6|), can be divided into three 
steps: (1) diffusive, (2) advective and (3) reactive. The 
first two steps determine the origin of the Lagrangian 
trajectory evolving with a given noise realization r) and 
ending in x. In the third step, the reaction at point x for 
the advected/diffused passive scalar 9 is computed: 



1) backward diffusion: x — » x — v / 2DAtrj 

2) backward advection by the Lagrangian map: 
x - y/2DAtrj -> F~ A *(x - y/2DAtri) 

3) forward reaction: 

6(t + At) = G At {6(t)). 

This three steps can be numerically implemented as fol- 
low. For each grid point x„ m , one uses N independent 
Gaussian processes W°, a = 1,...,N, N ^> 1, and 
computes x™ m = x„ >m — \/2D At W a . Then, using the 
Lagrangian backward propagator, r" m = F~ A *(x" m ). 
For 0r h m{t + At) one needs the values of 9 at time t in 
the positions r" m . Typically the r" are not on the 
grid points [nAx, mAy), nevertheless we can compute 
the value #(r" m , t) using simple linear interpolation from 
9n.m{t)- Therefore we have 

1 N 

6 n>m {t + At) = -Y,G[6{r« !m ,t)]. 



Similar to the previous case, one can integrate the dy- 
namics of the optical front using a two step discrete-time 
process. 
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FIG. 10. Pictorial scheme of the numerical algorithm for 
the geometrical optics limit. 



Starting from the field Q n ,m at time t, one can obtain the 
field at time t + At with the following algorithm 

1) using the Lagrangian propagator one evolves the 
interface between burned and unburned region; 

2) at each point of the evolved interface one constructs 
a circle of radius vo At, burning the points within 
the circles. 
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To numerically implement such an algorithm one can pro- 
ceed as follows: starting in a grid point, x„ jm , of the 
scalar field at time t + At one applies the backward evo- 
lution arriving at the point y = F~ A *x njm at the time 
t. In this point we construct the circle of radius vq At. 
If in this circle there is at least one burned point of the 
scalar field at time t, we fix (9(x„. m ; t + At) — 1 otherwise 
0(x n<m ;t + At) = 0. 

Also in this case we have to care about the radius of 
circle Vq At has to be much larger than the grid size Ax. 
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